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Abstract 

Unitary space-time modulation using multiple antennas promises reliable communi¬ 
cation at high transmission rates. The basic principles are well understood and certain 
criteria for designing good unitary constellations have been presented. 

There exist two important design criteria for unitary space time codes. In the 
situation where the signal to noise ratio is large it is well known that the diversity 
product of a constellation should be as large as possible. It is less known that the 
diversity sum is a very important design criterion for codes working in a low SNR 
environment. For some special situations, it will be more practical and reasonable to 
consider a constellation optimized at certain SNR interval. For this reason we introduce 
the diversity function as a general design criterion. So far, no general method to design 
good-performing constellation with large diversity for any number of transmit antennas 
and any transmission rate exists. 

In this paper we introduce numerical methods which allows one to construct codes 
with near optimal diversity sum. We also demonstrate how these numerical techniques 
lead to codes with excellent diversity product and more generally diversity function. 
The numerical design methods work for any dimensional constellation having any trans¬ 
mission rate. Codes can be optimized for any signal to noise ratio. 
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was also supported by a fellowship from the Center of Applied Mathematics at the University of Notre Dame. 
A preliminary version of this paper was presented at 40-th Allerton Conference on Communication, Control, 
and Computing, Monticello, Illinois, October 2002. 




1 Introduction and Model 


One way to acquire reliable transmission with high transmission rate on a wireless channel is 
to use multiple transmit or receive antennas. Either because of rapid changes in the channel 
parameters or because of limited system resources, it is reasonable to assume that both the 
transmitter and the receiver don’t know about the channel state information (CSI), i.e. the 
channel is non-coherent. 

In , Hochwald and Marzetta study unitary space-time modulation. Consider a wireless 
communication system with M transmit antennas and N receive antennas operating in a 
Rayleigh flat-fading channel. We assume time is discrete and at each time slot, signals are 
transmitted simultaneously from the M transmitter antennas. We can further assume that 
the wireless channel is quasi-static over a time block of length T. 

A signal constellation V := {$ 1 ,...,$^} consists of L matrices having size T x M and 
satisfying T > M and = Im- The last equation simply states that the columns of 

form a “unitary frame”, i.e. the column vectors all have unit length in the complex 
vector space and the vectors are pairwise orthogonal. The scaled matrices y/T^k, k = 
1, 2, • • • , L, represent the code words used during the transmission. One can verify that the 
transmission rate is determined by L and T: 


^ log2(T) 
T 


Denote by p the signal to noise ratio (SNR). The basic equation between the received 
signal R and the transmitted signal y/T^ is given through: 


R 



^H + W, 


where the M x N matrix H accounts for the multiplicative complex Gaussian fading coeffi¬ 
cients and the T x N matrix W accounts for the additive white Gaussian noise. The entries 
hm,n of the matrix H as well as the entries Wt^n of the matrix W are assumed to have a 
statistically independent normal distribution CA/'(0,1). In particular it is assumed that the 
receiver does not know the exact values of either the entries of if or W (other than their 
statistical distribution). 

The decoding task asks for the computation of the most likely sent code word given 
the received signal R. Denote by || \\f the Frobenius norm of a matrix. If A = (ajj) then 

the Frobenius norm is defined through 11A| Ij? = Under the assumption of above 

model the maximum likelihood (ML) decoder will have to compute: 


<hML = arg max ||R*<hz||p, 


for each received signal R. (See El 


Let (5m(‘h^<h;/) be the m-th singular value of It has been shown in El that the 

pairwise probability of mistaking for using Maximum Likelihood Decoding satisfies: 
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4(1 + pT/M) 


dtc (1.1) 
( 1 . 2 ) 


It is a basic design objective to construct constellations V = {<hi,..., such that the 
pairwise probabilities Plj.are as small as possible. Mathematically we are dealing with 
an optimization problem with unitary constraints: 


Minimize with the constraints = I where i = 1,2, 


,L. 


Formula (p..2|) is sometimes referred to as “Chernoff’s bound”. This formula is easy to 
work with, the exact formula O) is in general not easy to work with, although as we 
demonstrate in this paper it is very useful in the numerical search of good constellations 
as well. Researchers have been searching for constructions where the maximal pairwise 
probability of is as small as possible. Of course the pairwise probabilities depend on 

the chosen signal to noise ratio p and the construction of constellations has therefore to be 
optimized for particular values of the SNR. 

The design objective is slightly simplihed if one assumes that transmission operates at 
high signal to noise ratios. In []^, a design criterion for high SNR is presented and the 
problem has been converted to the design of a hnite set of unitary matrices whose pairwise 
diversity product is as large as possible. In this special situation several researchers ^ 
came up with algebraic constructions and we will say more about this in the next 

section. 

The main purpose of this paper is to develop numerical procedures which allow one to 
construct unitary constellations with excellent diversity for any set of parameters M, N, T, L 
and for any signal to noise ratio p. 

The paper is structured as follows. In Section || we introduce the diversity function of a 
constellation. This function depends on the signal to noise ratio and it gives for each value p 
an indication how good the constellation V will perform. For large values of p the diversity 
function is governed by the diversity product, for small values of p it is governed by the 
diversity sum. These concepts are introduced in Section || as well. The introduced concepts 
are illustrated on some well known constellations previously studied in the literature. 

Section ^ is probably the main section of this paper. In this section we demonstrate 
how one can use numerical tools to derive excellent constellations for any set of parameters 
M, N, T, L and p. The main numerical techniques we are using are the Simulated Annealing 
Algorithm and the Genetic Algorithm which are both explained in Section 

In Subsection H we present a larger set of 2-dimensional code designs and we indi¬ 
cate their performances through the computation of their diversity and through simulations. 
Subsection gives results for general n-dimensional constellations. The simulations indi- 
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cate that in the design of codes more attention should be given to the diversity sum (more 
generally diversity function) which previously has not been studied. 


2 The Diversity Function, the Diversity Sum and the 
Diversity Product 


In this paper we will be concerned with the numerical construction of constellations where 
the right hand sides in (|1 . 1|) and ( p..2|) , maximized over all pairs 1,1' is as small as possible 
for hxed numbers of T, M, N, L. As already mentioned this tasks depends on the signal to 
noise ratio the system is operating. For this purpose we define the exact diversity function 
dependent on the constellation V = {<Fi,..., $^1 and a particular SNR p through: 


Veiy,p) 


max Prob ( choose <F;/ 


transmitted ) (p) 


( 2 . 1 ) 


For a particular constellation with a large number L of elements, with many transmit 
and receive antennas the function T>e{V,p) is very difficult to compute. Indeed for each pair 
<!>;/, it is required to compute the singular values of the M x M matrix and then 

one has to evaluate up to L{L — l)/2 integrals of the form (|l.lj) and this has to be done for 
each value of p. Although this task is formidable it can be done in cases where T, M, L are 
all in the single digits using e.g. Maple. 

Using Chernoff’s bound (|1.2|) we dehne a simplified function called the diversity function 
through: 


M r 


Viy,p) := max- 


n 

m=l 


1 + 


{pTiMf 
4(1 + pT/M) 




-N 


( 2 . 2 ) 


The computation of V{V,p) does not require the evaluation of an integral and the com¬ 
putation requires essentially the computation of ML{L — l)/2 singular values. 

The singular values 5m(‘h^4);/) are by dehnition all real numbers in the interval [0,1] 
as we assume that the columns of form both orthonormal frames. The functions 

Ve{y,p) and V{y,p) are smallest if the singular values are as small as possible. 

These numbers are all equal to zero if and only if the column spaces of <F;, are pairwise 
perpendicular. We call such a constellation fully isotropic. Since the columns of generate 
a M-dimensional subspace this can only happen if L < T/M. On the other hand ii L <T/M 
it is easy to construct a constellation where the singular values of (<I>^<F;/) are all zero. Just 
pick LM different columns from a T x T unitary matrix. Figure ^ depicts the functions 
Ve{y,p) and V{y,p) for a fully isotropic constellation with T = 10 and M = N = 2. 

In order to study the function T>{V,p) more carefully let 


(pT/M) 
i{l+pT/M)' 


In some small interval [pi, p 2 \ the maximum in (p.2|) is achieved for some hxed indices 1,1' 
and in terms of p the function TPiV, p) is of the form: 


v{y,p) 


1 

2 (1 -f- cipH-h cmP^)^' 
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Figure 1: Diversity function D(V, p) and exact diversity function De(V, p) of a fully isotropic 
constellation. 

where the coefficients Ci,... ,cm depend on the particular constellation and on the chosen 
interval [pi, p^. For an interval close to zero the dominating term will be the coefficient Ci. Up 
to some factor this term will dehne the diversity sum of the constellation. When p >> 0 then 
the dominating term will be the coefficient cm and up to some scaling this term will dehne 
the diversity product of the constellation. A constellation will have a small diversity function 
for small values of p (and presumably performs well in this range) when the constellation 
is chosen having a large diversity sum. A constellation will have a small diversity function 
for large values of p (and presumably performs well in this range) when the constellation 
is chosen having a large diversity product. In the next two subsections we will study the 
limiting behavior of T>{V,p) as p goes to zero and to inhnity. 

2.1 Design criterion for high SNR 

When the SNR p is very large then D(V, p) can be approximated via: 

- 1*1' 2 V4(l + pT/M)) n (1 ^ 

It is the design objective to construct a constellation <I)i, $ 2 , • • ■ , such that 

M 
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is as large as possible. This last expression defines in essence the diversity product. In order 
to compare different dimensional constellations it is customary to use the definition: 

Definition 2.1. (See |]^) The diversity product of a unitary constellation V is defined as 

( M 

m=l 

An important special case occurs when T = 2M. In this situation it is customary to 
represent all unitary matrices in the form: 



^k = 


2 



(2.5) 


Note that by definition of the matrix is a M x M unitary matrix. The diversity 
product as defined in Definition p.l| has then a nice form in terms of the unitary matrices. 
For this let Am be the mth eigenvalue of a matrix. Then we have: 


1 - 


— ^Am(2/7vr 


^*^1, 






M 








So we have 


M 

m=l 


^ n ~ = ^1 det(Tz. - 

m=l 


When T = 2M and the constellation V is defined as above, then the formula of the diversity 
product assumes the simple form: 



- min I det(\I'i — . 

2o<i<i'<L 


( 2 . 6 ) 


We call a constellation V a fully diverse constellation if V > 0. 

A lot of efforts have been taken to construct constellations with large diversity product. 
(See e.g. |^, |T^, |^, For the particular situation T = 2M the design asks for 

the construction of a discrete subset V = {\ki,...,T 2 ,}of the set of M x M unitary matrices 
U{M). When this discrete subset has the structure of a discrete subgroup of U{M) then the 
condition that V is fully diverse is equivalent with the condition that the identity matrix is 
the only element of V having an eigenvalue of 1. In other words the constellation V is required 
to operate fixed point free on the vector space C^. Using a classical classification result of 
fixed point free unitary representations by Zassenhaus Shokrollahi et al. [^, were 


able to study the complete list of fully diverse finite group constellations inside the unitary 
group U{M). Some of these constellations have the best known diversity product for given 
fixed parameters M, N, L. Unfortunately the possible configurations derived in this way is 
somehow limited. The constellations are also optimized for the diversity product and as we 
demonstrate in this paper equal attention should be given to the diversity sum. 
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Remark 2.2. In most of the literature mentioned above researchers focus their attention to 
constellations having the special form (p.5|) . Unitary differential modulation is used to avoid 
sending the identity (upper part of every element in the constellation) redundantly. This 
increases the transmission rate by a factor of 2 to: 

log2(-^) ^ o log2(-^) 

M T 

Because of this reason we will also focus ourselves in the later part of the paper to the special 
form (p.5|) as well. Nonetheless it will become obvious that the numerical techniques also 
work in the general situation. 


2.2 Design criterion for low SNR channel 


As we mentioned before a constellation with a large diversity sum will have a small diversity 
function at small values of the signal to noise ratio. This is particularly suitable when the 
system operates in a very noisy environment. When p is small, using Formula ( ^.3| ), one has 
the following expansion: 


n =n ii+ 

m=l \ ^ r I ) 

M 

= (1 + p)“ - p ^ + p)"-‘ + 0(if). 

m=l 

When p —0, i.e. p —0, we can omit the higher order terms 0{p^) and the upper bound 
of requires that 

m 

112 

is small. In order to lower the pairwise error probability, it is the objective to make || <!>)'$;/||p, 
as small as possible for every pair of /, V. It follows that at high SNR, the probability primarily 
depends on nm=i(l — S'^), but at low SNR, the probability primarily depends on J2m=i 
In order to be able to compare the constellation of different dimensions, we dehne: 

Definition 2.3. The diversity sum of a unitary constellation V is dehned as 



min 


f- 


Again one has the important special case where T 
special form (|2.5|) . In this case one verihes that 


2M and the matrices take the 


l^z = ^11^ + ^z = ^tr ((/ + T;,Tz)(J + 

= ^tr {21 + = ^(4M - (2M - tr 


= i(4M - tr ((T^ - T,)*(^z - ^zO)) = 


IT, - T,/ 


i'Wf) 
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For the form ( p.5|) the diversity sum assumes the following simple form: 

1 




= mm 


t;' 2yM 




I'WF- 


(2.7) 


Remark 2.4. Without mentioning the term the concept of diversity sum was used in [^. 
Liang and Xia |^, p. 2295] explicitly defined the diversity sum in the situation when 
T = 2M using equation (|2.7|) . Dehnition |2.3| naturally generalizes the dehnition to arbitrary 
constellations. 

As the formulas make it clear the diversity sum and the diversity product are in general 
very different. There is however an exception. When T = 4, M = 2 and the constellation V 


is in the special (^). If in addition all the 2x2 matrices {Ti,..., are a subset of the 
special unitary group 


SU{2) = {Ae \a*A = I and det A = 1} 

then it turns out that the diversity product ]]][ V and the diversity sum ^ V of such a 
constellation are the same. For this note that elements of SU{2) have the special 

form: 



Through a direct calculation one verihes that det(4/i — T;/) = |a — cp + |6 — dp and 
||4/; — = 2(|a — cp + \b — dp). But this means that ]Q V = for constellations 

inside SU{2). 

2.3 Three illustrative examples 

The diversity sum and the diversity product govern the diversity function at low SNR re¬ 
spectively at high SNR. Codes optimized at these extreme values of the SNR-axis do not 
necessarily perform very well on the “other side of the spectrum”. In this subsection we illus¬ 
trate the introduced concepts on three examples. All examples have about equal parameters, 
namely T = 4, M = 2 and the size L is 121 respectively 120. The hrst two examples are 
well studied examples from the literature. The third example is one we derived by numerical 
methods. 

Orthogonal Design: This constellation has been considered by several authors p0| . 
For our purpose we simply define this code as a subset of SU{2): 

2rmri 2mTi \ 

2 11 e 11 \ I 

2mfi _ 2m7ri I |'m,,'n. = 0,l,-- - ,10 

e 11 e 11 / 

The constellation has 121 elements and the diversity sum and the diversity product are the 
same. 
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Unitary Representation of S'L 2 (F 5 ): Shokrollahi e. a. derived a constellation using 
the theory of hxed point free representations whose diversity product is near optimal. This 
constellation appears as a unitary representation of the hnite group S'L 2 (F 5 ) and we will 
refer to this constellation as the S'L 2 (F 5 )-constellation. The hnite group S'L 2 (F 5 ) has 120 
elements and this is also the size of the constellation. In order to describe the constellation 
let rj = e~ and dehne 


P = 


2 2 1 
T] — rj rj 


rj 


^\rj^ -rj^ rj^ -Tj'^ 


, Q = 


1 2 9 1 

Tj — rj Tj — Tj 


a /5 \ rj^ — Tj^ rj^ — rj^ 


Then the constellation is given by the set of matrices {PQyX,where j = 0,1, ■ ■ ■ ,9, X runs 
over the set 


{h, P, Q, QP, QPQ, QPQP, QPQ\ QPQPQ, QPQPQ\ 

QPQPQ^P, QPQPQ^PQ, QPQPQ^PQP}. 

The constellation has rate R = 3.45 and H 5'T2(F5) = ^ S'L 2 (F 5 ) = ~ 0.309. 

The diversity product of this constellation is truly outstanding. For illustrative purposes we 
plotted in Figure ^ the diversity functions of this constellation. 
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Figure 2: Diversity function D(V,p) and exact diversity function for the group constellation 
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Numerically Derived Constellation: Using Simulated Annealing Algorithm we found 
after short computation a constellation with near perfect diversity sum. The constellation 
is given through a set of 121 matrices 




k,l 


B = 


/ -0.9049+ 0.3265 0.1635 + 0.2188 * z \ 

V 0.0364 +0.2707 *z -0.8748 + 0.4002 * i 


-0.1596 + 0.9767 * i -0.1038 + 0.0994 * i 
0.0833 -0.1171 *z -0.9432 +0.2995 *z 





Table 1. The following table summarizes the parameters of the three constellations: 



Orthogonal 

design 

5L2(F5) 

Numerically 

derived 

Number of elements 

121 

120 

121 

diversity sum 

0.1992 

0.309 

0.3886 

diversity product 

0.1992 

0.309 

0.0278 


Of course we were curious about the performance of these three different codes. Figure ^ 
provides simulation results for each of the three constellations: 



Figure 3: Simulations of three constellations having sizes T = 4, M = 2 and L = 120 
respectively L = 121. 
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Note that the numerically designed code who has a very bad diversity product is per¬ 
forming very well nevertheless due to the exceptional diversity sum. These simulation results 
give an indication that the diversity sum is a very important parameter for a constellation. 


3 Numerical Design of Unitary Constellations with Good 
Diversity 

In order to numerically design constellations with good diversity it will be necessary to 
have a good parameterization for the set of all constellations having size L, block length 
T and operating with M transmit antennas. In this section we show how one can use the 
theory of complex Stiefel manifolds and the classical Cayley transform to obtain such a 
parameterization in all cases. 


3.1 The complex Stiefel manifold 

Definition 3.1. The subset of T x M complex matrices 

St,m ■= e I = Im} 

is called the complex Stiefel manifold. 

From an abstract point of view a constellation V := {<Fi,..., <hL} having size L, block 
length T and operating with M antennas can be viewed as a point in the complex manifold 

M. := {St,m)^ = St,m X • • • X St,m ■ 

L copies 

The search for good constellations V requires hence the search for points in At whose diversity 
is excellent in some interval [pi,p 2 ]- 

Stiefel manifolds have been intensly studied in the mathematics literature since their 
introduction by Eduard Stiefel some 50 years ago. A classical paper on complex Stiefel 
manifolds is 0, a paper with a point of view towards numerical algorithms is ^j. The major 
properties are summarized by the following theorem: 

Theorem 3.2. St, M is a smooth, real and compaet sub-manifold of of real 

dimension 2TM — M^. 


Some of the stated properties will follow from our further development. 
The following two examples give some special cases. 

Example 3.3. 


St,i = U ^ I 11^11 ~ 

\ ~ 

is isomorphic to the 2T — 1 dimensional unit sphere 


M 


^ ^ XiXi = 1 ^ c 


t,2T 
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Example 3.4. When T = M then St,m = U{M), the group of M x M unitary matrices. It 
is well known that the Lie algebra of U{M), i.e. the tangent space at the identity element, 
consists of all M X M skew-Hermitian matrices. This linear vector space has real dimension 
in particular the dimension of U{M) is as well. 

A direct consequence of Theorem |3.2| is: 

Corollary 3.5. The manifold M which parameterizes the set of all constellations V having 
size L, block length T and operating with M antennas forms a a real compact manifold of 
dimension 2LTM — . 


As this corollary makes it clear a full search over the total parameter space is only possible 
for very moderate sizes of M, L, T. It is also required to have a good parameterization of 
the complex Stiefel manifold St,m and we will go after this task next. 

The unitary group is closely related to the complex Stiefel manifold and the problem of 
parameterization ultimately boils down to the parameterization of unitary matrices. For this 
assume that <F is a T x M matrix representing an element of the complex Stiefel manifold 
St,m- Using Gramm-Schmidt one constructs aT x [T — M) matrix V such that the T xT 
matrix [*F | U] is unitary. Dehne two T x T unitary matrices [<Fi | Vi] and [<F 2 | V 2 ] to be 
equivalent whenever $1 = 4 ) 2 . A direct calculation shows that two matrices are equivalent 
if and only if there is (T — M) x (T — M) matrix Q such that: 


[<h2 I U 2 ] 




(3.1) 


Identifying the set of matrices Q appearing in (|3.1|) with the unitary group U{T — M) we 
get the result: 


Lemma 3.6. The complex Stiefel manifold St,m is isomorphic to the guotient group 


U{T)IU{T - M). 


This lemma let us verify the dimension formula for St,m 


stated in Theorem 3.2 


dim 5 t,m = dimG(r) - dimG(r - M) = T^ - (T - Mf = 2TM - . 


The section makes it clear that a good parameterization of the set of constellations 
V requires a good parameterization of the manifold M. and this in turn requires a good 
parameterization of the unitary group U{M). 

Once one has a nice parameterization of the unitary group U (M) then Lemma |3.6| pro¬ 
vides a way to parameterize the Stiefel manifold St,m as well. Parameterizing U{T) modulo 
U{T — M) is however an ‘over parameterization’. Edelman, Arias and Smith P] explained a 
way on how to describe a local neighborhood of a (real) Stiefel manifold St,m- The method 
can equally well be applied in the complex case. We do not pursue this parameterization in 
this paper and leave this for future work. 

In the remainder of this paper we will concentrate on constellations having the special 
form (p.5|) . From a numerical point of view we require for this a good parameterization of 
the unitary group and the next subsection provides an elegant way to do this. 


12 




3.2 Cayley transformation 


There are several ways to represent a unitary matrix in a very explicit way. One elegant way 
makes use of the classical Cayley transformation. In order that the paper is self contained 
we provide a short summary. More details are given in [^, Section 22] and ||^. 


Definition 3.7. For a complex M x M matrix Y which has no eigenvalues at —1, the Cayley 
transform of Y is dehned to be 


Y‘^={I + Y)-\I-Y), 
where I is the M x M identity matrix 

Note that (/ + Y) is nonsingular whenever Y has no eigenvalue at -1. One immediately 
verihes that {Y'^Y = Y. This is in analogy to the fact that the linear fractional transformation 
/(^) = the property that f{f{z)) = z. 

Recall that a matrix M is skew-Hermitian whenever A* = —A. The set oi M x M 
skew-Hermitian matrices forms a linear subspace of having real dimension 

M^. This is the Lie algebra of the unitary group U{M). The main property of the Cayley 
transformation is summarized in the following theorem. (See e.g. [^, p!B|). 

Theorem 3.8. When A is a skew-Hermitian matrix then (/ -|- A) is nonsingular and the 
Cayley transform V ■= A^ is a unitary matrix. Vice versa when V is a unitary matrix which 
has no eigenvalues at —1 then the Cayley transform is skew-Hermitian. 

This theorem allows one to parameterize the open set of U{M) consisting of all unitary 
matrices whose eigenvalues do not include —1 through the linear vector space of skew- 
Hermitian matrices. 

The Cayley transformation is very important for the numerical design of constellations 
because it makes the local topology of U{M) clear. One can see that most optimization 
method require us to consider the neighborhood of one element in U{M). 

3.3 Simulated Annealing Algorithm (SA) 

In our numerical experiments we have considered several methods. Because there is a large 
number of target functions the best known optimization algorithms such as Newton’s Meth¬ 
ods [T^, 1^ and the Conjugate Gradient Method are difficult to implement. 

Two algorithms of very general nature, the Simulated Annealing Algorithm and the Ge¬ 
netic Algorithm turned out to be very practical. In this and the next subsection we describe 
these algorithms. 

Simulated Annealing is a method which mimics the process of melted metal getting cooled 
off. In the annealing process of the melted metal, first the metal is heated to melt, then the 
temperature is getting down gradually. The metal will get to a minimized energy state if 
the temperature is lowering slow enough. For more details about this algorithm, we refer 
to [|, 

In fact, we would rather call it a general method instead of a concrete algorithm. Gen¬ 
erally speaking, for a given optimization problem, we always take an initial solution in some 
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certain way. Then consider the “neighborhood” of this solution, we will accept the solu¬ 
tion in the “neighborhood” according to some predefined criterion which might include a 
probability threshold. 

SA is a stochastic method which will hnd better solution up to the optimal point in 
a iterative way. And with Cayley transformation which is a good representation of any 
dimensional unitary matrix, one can see that this method can be applied to any dimensional 
constellation design based on the algebraic structure we proposed as above. 

The algorithm we are using also depends on the initial guess. In 2 dimensional constel¬ 
lation search, one can use the Brute Force with coarser grid to find a good guess, then use 
SA to iteratively find better and better points. But when it comes to higher dimension, it is 
not feasible to use this approach anymore with limited computing power. One way to find a 
good initial guess is to start from an existing constellation (such as a group constellation). 

Our implementation of the algorithm can be summarized in the following way: 

1. Generate initial generators of the whole constellation. For low dimension, one can try 
to use coarser mesh to hnd better starting points, for high dimension one can randomly 
choose an initial guess. 

2. Using the Cayley transformation of the generators generate randomly a new constel¬ 
lation in the neighborhood of the old constellation where the selection is done using a 
Gaussian distribution with decreasing variances as the algorithm progresses. 

3. Calculate the diversity function (product, sum) of the newly constructed constellation. 

4. If the new constellation has better diversity function (product, sum), then accept the 
new constellation. If not, reject the new constellation and keep the old constellation 
(or accept it according to Metropolis’s criterion [|T^). 

5. Check the stopping criterion, if satished, then stop, otherwise go to 2 and continue the 
iteration. 


3.4 Genetic Algorithm (GA) 


Genetic Algorithm [Q is a optimization algorithm proposed by J.H. Holland emulating 
the evolutionary process of species. GA doesn’t assume much specific information about 
the given problem. First a group of candidates are selected in a certain way and they are 
encoded using binary coding in most cases. Then consider the offsprings of the candidates 
(the mating processing is also defined with respect to the specific problem). In the same time 
there might be new random candidates added in. The whole idea is “better survive”, which 
really means some candidates with higher Fitness Evaluation Function keep staying in the 
group and other “worse” ones will get discarded. The algorithm will stop if some threshold 
is reached after certain number of iterations. The interested reader will find further details 


m 


Our adapted algorithm is summarized in the following way: 


1. Generate the initial population randomly with the desired size. 
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2. For every individual in the population, calculate its Fitness Evaluation Function. In our 
situation this is the minimum distance from this individual to all other individuals. Of 
course this depends on the criterion we want to optimize (diversity function, product, 
sum). 

3. Replace the individuals which have the worst Fitness Evaluation Function with the 
same number of randomly chosen individuals. (If necessary, randomly choose certain 
number of individuals for mutation). 

4. Calculate the diversity function (product, sum) of the given size population, if it im¬ 
proves, accept it, otherwise, restore the population to the previous population. 

5. Check the stopping criterion, if satished, then stop, otherwise go to 2 and continue the 
iteration. 

3.5 Constellations with algebraic structure 

As explained in the beginning of this section the number of parameters to be optimized 
explodes for even moderately chosen numbers T, M and L. The full parameterization of a 
constellation with size L = 120, T = 4 and M = 2 has 1440 free parameters and this is out 
of the range for a naive implementation of the simulated annealing algorithm. The genetic 
algorithm is still feasible in this range. 

In this subsection we explain how one can restrict the parameter space to judiciously 
chosen subsets and do the optimization inside this restricted parameter space only. 

Consider a general constellation of square unitary matrices: 

V = {Ti,T2,-- - ,T4. 

In order to calculate the diversity product (or sum), one needs to do calculations: 

I det(\I'j — 'hj)|, for every different pair i,j. 

If one deals with a group constellation then one needs only to calculate L — 1 such 
determinant calculations and this is one of the big advantages of group constellations. This 
is a direct consequence of 

I det(Ti - Tj)| = I det(Ti) det(/ - T*T^)| = | det(/ - 
where is still in the group. 

Group constellations are however very restrictive what the algebraic structure is con¬ 
cerned. In the following we are going to present some constellations in between general 
constellations and group constellations. Our constellations have some small number of gen¬ 
erators. This will ensure that the total parameter space to be searched is limited as well. 
Like for group constellations we can reduce the number of targets. We start with an example: 

Example 3.9. Consider the constellation 

V = {A’‘B^\A,B eU{M),k = 0,--- ,p,/ = 0,--- ,g}. 
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The parameter space for this constellation is U{M) x U{M), this is a manifold of dimen¬ 
sion 2 M^ and the number of elements is {p+ !)((?+ 1). If one has to compute | det(Tj — Tj)! 


for every distinct pair this would require 


We will show in the following 


(p+ l)(g + 1) 

2 

that the same result can be obtained by doing 2pq + p + q determinant computations. This 
is in analogy to the situation of group constellations. 

Let Tj and Tj be two distinct elements having the form and respectively. 

We have now several cases. When ki = k 2 , then necessarily li ^ I 2 and the distance is 
computed as 

I det(v4^i5'i - A'^^B^^)\ = I det(/ - 

where {h — h\ is an integer between 1 and q. 

If /i = I 2 , then we have ki 7 ^ k 2 and the distance is computed as 


det{A’^^B^^ - A'^^B^^)] = I det(J - 


where \k2 — ki\ is an integer between 1 and p. 

If (fci < ^2 and li < I2) or {ki > ^2 and li > I2), we have 

I det{A^^B^^ - A’^^B^^) \ = \ det(/ - A^k2-k^\^\i2-h\^^ 

where 1 < 1^2 — fci| < p and 1 < {h — h\ < <?• 

Similarly if (fci < ^2 and h > I2) or {ki > k2 and h < I2) then 

I det(/l^iS'i - = I det(Lll^^-^il - 


with 1 <\k 2 — ki\<p and 1 < \l 2 — li\ < P- 

The total number of distances to be computed is in total equal to 2pq + p + q. 

Above idea can be applied to describe in an algebraic way constellations where the number 
of distances is always considerably smaller than the total number of comparisons involving 
all possible pairs. The following provides a list of similar generator sets for constellations. 


1. V = {A^B^\A,B eU{M),k = d,--- ,p,/= O,-- - ,g} 

2. V = B,C E U{M), k = 0, - ■ ■ ,p,l = 0, ■ ■ ■ ,q,m = 0, - ■ ■ , r} 

3. I,A,AB,ABA,ABAB,ABABA,--- 

4. /, A, AB, ABC, ABC A, ABCAB, ABC ABC, ABC ABC A, ■ ■ ■ 

5. I,AB,A‘^B\A^B^,--- 

6 . I,ABC,A^B^C^,A^B^C^,--- 


The last generator set can also be found in |13 


to design codes algebraically . Above 
construction ideas can be applied to do “product constructions”. For instance, 
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7. Let Si = {/, (7, (7^, C^, • • ■ } and S 2 = {I, A, AB, ABA, • • • } and consider the cartesian 
product constellation 


S — Si X S 2 — {SiS2|Si G 5*1,52 G 5’2}. 


8. Similarly, let = {/, A, AB, ABA, • ■ ■ } and ^2 = {/, C, CD, DCD, CDCD, • ■ ■} and 
consider the constellation 

S' = 5*1 X 5*2 = {SiS2|Sl G Si,S2 £ 5'2}. 

For a P, P' G U{M) and any matrix Q G , we have 

||PQP'il^ = ||g||^. 

Like for the diversity product it follows from above that the number of targets to be checked 
in order to compute the diversity sum is reduced in a similar way. In fact, it can be shown 
that the method works in the same way if one wants to optimize the diversity function at a 
certain SNR. This is essential the case because: 

SUPQn = SmiQ), ( 3 . 2 ) 


for m = 1, 2, • • • , M. 

With constellations as above we are able to reduce the dimension of the parameter space 
and in the same time we have a considerable reduction in the number of targets to be checked. 
This is absolutely essential for the Simulated Annealing Algorithm which we described before. 


Remark 3.10. Because of (|3.2|) one can use the algebraic structures we presented also for 
the design of constellations where exact diversity function is optimized at a certain SNR. 
The idea is appealing but of course it involves the evaluation of an improper integral in every 
loop and this is computationally very expensive. For a constellation of size 3 this is however 
easily possible and we have done this in the sequel. 

The hrst constellation is 3 element 2 dimensional constellation with optimal diversity 
product and diversity sum(consult Appendix ^ for more details): 


I,A,B\A = 


‘I'Ki 

e 3 


— 27r't 

e 3 


,5 = 



The second constellation is optimized at 5 db by Simulated Annealing Algorithm based 
on exact diversity function: 


/ -0.4530000 - 0.7804689 * i, 0.2197119 - 0.3706563 * i \ 

V -0.1733377 - 0.3944787 * i, -0.5439448 + 0.7200449 * i j ’ 


/ -0.4475155 + 0.7358245 * i, -0.1243078 + 0.4927877 *i\\ 
V 0.1862253 +0.4728766 -0.5378253 - 0.6726454 * z Jj’ 


The simulation results are given in Figure ^ 
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Figure 4: Optimized constellation by exact diversity function versus constellation with op¬ 
timal diversity product and sum. From the graph, one can see that the two constellation 
almost have the same performance, except there are minor deviations at certain interval. 
Actually it can be verihed that the elements in the second constellation are similar to (up 
to one unitary matrix) those in the first constellation, so it has also optimal diversity prod¬ 
uct and diversity sum. For more details about 2 dimensional optimal constellation with 3 
elements, one can look at Appendix 


4 Numerical Results of Constellation Design 

4.1 2 Dimensional constellation design 

We are going to talk about 2 dimensional constellation design separately because unlike 
higher dimensional ones, we have a very explicit representation of U{2): 

^ ||a |2 +| 6|2 = 1 , O< 0 < 27 r|. (4.1) 

From above, one can see that U{2) is 4 dimensional, something we have already mentioned 
before. So when we consider the hnitely generated constellation with a good algebraic 
structure, it is possible to use Brute Force method to calculate the almost optimal generator. 
Of course in this process, hner grid and more computing power will make better result. 


U(2) = 


ae* 


-b 
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Table 2. Comparisons of different methods and different parameterizations: 



size 

diversity product 

structure 

Simulated Annealing 

125 

0.2127 

~WWc^ 

Simulated Annealing 

120 

0.2202 

WW 

Simulated Annealing 

121 

0.2417 


Brute Force 

120 

0.1914 


Genetic Algorithm 

120 

0.2377 

N/A 


Table 3. Optimization of the the diversity product using Simulated Annealing and different 
generator sets. All computations were accomplished in less than 3 minutes. 


structure/size 

36 

49 

64 

256 

400 

900 

10000 


0.3860 

0.3781 

0.2742 

0.1025 

0.0866 

0.0834 

0.0158 

AB 

0.3205 

0.2659 

0.2450 

0.1030 

0.0800 

0.0579 

0.0122 

7^ 

0.3769 

0.3502 

0.3090 

0.1651 

0.1342 

0.0820 

0.0187 


structure/size 

27 

64 

216 

343 

512 

729 

9261 

WWc^ 

0.3418 

0.2616 

0.1833 

0.1401 

0.0632 

0.1012 

0.0031 

ABC 

0.3299 

0.1832 

0.1033 

0.0725 

0.0555 

0.0430 

N/A 


0.4122 

0.2512 

0.0583 

0.0206 

0.0087 

0.0039 

N/A 


The above tables show some of the numerical results in 2 dimension. The constellations 
in Table 1^ all have about 120 elements and they are therefore comparable with the S'L 2 (F 5 ) 
constellation described in Subsection ^7^ . The obtained diversity products are reasonable 
good but fall short from the outstanding diversity product of 0.309 which the S'L 2 (F 5 ) 
constellation has. 

The numerical design method seems to be particular powerful if one seeks constellations in 
higher dimension or if one wants to optimize the diversity sum. One can see the performance 
of higher dimensional constellations in Subsection [4.3i 

Table 4. Comparison of different methods or different parameters 



size 

diversity sum 

structure 

Simulated Annealing 

125 

0.3919 

~WWc^ 

Simulated Annealing 

120 

0.3696 

WW 

Simulated Annealing 

121 

0.3886 

WW 

Brute Force 

120 

0.3673 

WW 

Genetic Algorithm 

120 

0.3867 

N/A 
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Table 5. Optimization of the the diversity sum using Simulated Annealing and different 
generator sets. All computations were accomplished in less than 3 minutes. 


structure/size 

36 

49 

64 

256 

400 

900 

10000 


0.5113 

0.4733 

0.4474 

0.2875 

0.2504 

0.1848 

0.0785 

AB 

0.5530 

0.4240 

0.3821 

0.1994 

0.1629 

0.1064 

0.0310 


0.5466 

0.5121 

0.4735 

0.3088 

0.2637 

0.2047 

0.0869 


structure/size 

27 

64 

216 

343 

512 

729 

9261 

WWc^ 

0.5400 

0.4210 

0.2992 

0.2663 

0.2099 

0.2060 

0.0772 

ABC 

0.5382 

0.4497 

0.2614 

0.2065 

0.1695 

0.1447 

0.0398 


0.5630 

0.4271 

0.2864 

0.2198 

0.1969 

0.1423 

N/A 


Numerical method is extremely good for designing constellation with large diversity sum, 
most of the results are the best or very close to the best codes. It can be verihed that the 
S'L 2 (F 5 ) group code also has diversity sum 0.309, we believe that would be the reason its 
performance is so good. Again as for higher dimension and large cardinality, numerical 
methods are even more powerful. 

We have tried SA on constellation with 10000 elements using the structure A^BK For 
this we can get a diversity sum of 0.1000. For the construction in [[^] which can be proven to 
be asymptotically the best as a subset of SU{2), we can only get a diversity sum of 0.0654 
with 8433 elements or a diversity sum of 0.0604 with 10770 elements. 


4.2 Constellation design for any dimension 


As we mentioned before, we can always choose an existing constellation as staring point as 
our numerical method. In the sequel, we use the group constellation ^ 21,4 in p0| : 


f 7] 0 0 \ /010\ 

Vi = {A’^B^\A= \ 0 rj^ 0 ,5= 0 0 1 , A; = 0,1, • ■ ■ , 20, / = 0,1, 2} 

\ 0 0 77^^ / \ r ]'^ 0 0 J 

One can verify that 

JJVi = 0.3851 

It seems like the ^21,4 is already very good constellation, our algorithm only improves a 
little. However one can check most of the case, the algorithm will improve much compared 
to the original group constellation. 

V2 = {A^B^\k = 0 , 1 ,- ■■ ,20,/= 0,1, 2}, 

where 

/ 0.9415 + 0.3155 * i 0.0573 - 0.0222 * i 0.0496 + 0.0882 * i \ 

A=i 0.0160 - 0.0555*7 0.4005 + 0.9136*7 0.0326 - 0.0212 *7 , 

\ 0.0579 + 0.0855*7 -0.0312 -0.0099*7 0.1384-0.9844*7 / 
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B 


0.0175 +0.0095 
0.0086 +0.0100 
-0.4836 + 0.8750 * i 


0.9997 +0.0111 
-0.0082 + 0.0040 * i 
0.0004- 0.0198 


0.0079 + 0.0042 * i \ 
0.9999 + 0.0036 * i 
-0.0045 - 0.0126 / 


One verifies that 


JJVs = 0.3874. 

Another approach is to choose the starting points randomly. We have the following tables 
to show the results of some experiments. The experiments are based on the computation on 
a Intel Pentium 800MHz PC and no computation lasts longer than 3 minutes. 


Table 6. Numerical results on diversity product using Simulated Annealing Algorithm 


number of 
elements 

dim=2 

dim=3 

dim=4 

dim=5 

4 

0.7071 

0.7657 

0.7388 

0.6768 

9 

0.5701 

0.5754 

0.4774 

0.4259 

16 

0.4018 

0.4574 

0.4651 

0.3877 

25 

0.3443 

0.3834 

0.3809 

0.3467 

36 

0.2865 

0.3450 

0.3501 

0.3760 


Table 7. Numerical results on diversity sum using Simulated Annealing Algorithm 


number of 
elements 

dim=2 

dim=3 

dim=4 

dim=5 

4 

0.8147 

0.8160 

0.7861 

0.7377 

9 

0.6956 

0.6861 

0.6539 

0.6389 

16 

0.5908 

0.6459 

0.6288 

0.5916 

25 

0.5618 

0.6268 

0.6190 

0.5795 

36 

0.5286 

0.6054 

0.6148 

0.5853 


From above table, SA works well when the cardinality of the constellation is small. 
When it comes to large number of constellation, SA works even better as we mentioned 
above. Also SA’s results on higher dimension showing the pretty much the same size as 2 
dimension which makes us believe it will also work in higher dimension, although we don’t 
have many comparisons available. 

We have the following tables to show the results of some experiments using Genetic 
Algorithm. The experiments are based on the computation on a Intel Pentium 800MHz PC 
and no computation lasts longer than 3 minutes. 

Table 8. Numerical results on diversity product using Genetic Algorithm 


number of 
elements 

dim=2 

optimal DP 

dim=3 

dim=4 

dim=5 

3 

0.8644 

^=0.8860 

0.8264 

0.7305 

0.6737 

4 

0.8051 

y^=0.8165 

0.7343 

0.6521 

0.6305 

6 

0.6924 

unknown 

0.6632 

0.6154 

0.5721 

10 

0.5768 

unknown 

0.5497 

0.5742 

0.4942 
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Table 9. Numerical results on diversity sum using Genetic Algorithm 


number of 
elements 

dim=2 

optimal DS 

dim=3 

dim=4 

dim=5 

3 

0.8601 

^=0.8860 

0.8331 

0.8118 

0.7798 

4 

0.8029 

^1=0.8165 

0.7802 

0.7757 

0.7492 

6 

0.7443 

./|=0.7746 

V ^ 

0.7502 

0.7293 

0.7176 

10 

0.6826 

^=(0.7071) 

0.6981 

0.6920 

0.6817 


Genetic algorithm works extremely well when it comes to small number constellation. 
For instance, look at 3 elements 2 dimension diversity product case, one can check the 
optimal diversity product is ^ ~ 0.8660(see Appendix while our results shows 0.8601. 
The results make us believe for higher dimension, maybe we are approaching the optimal 
diversity product (sum) too (In the above table for 3 or higher dimension constellation, we 
don’t know the optimal diversity product). One can see other comparisons with known 
optimal diversity product (sum) |T^. Whereas when it comes to huge size constellation, 
our current GA doesn’t seems to work as well as small number constellation because of 
drastically increased complexity. Maybe some more complicated mechanism is needed to 
add in to handle this problem more subtly. 

4.3 Performance of different dimension constellation 

Basically our numerical method can be applied to design constellation for every dimension 
and any transmission rate. The following graph shows the comparison of three constellations 
with different dimensions with 2 receiver antennas. The hrst one is a 2 dimensional con¬ 
stellation with 3 elements {R = 0.7925) and optimal diversity product 0.8660 and optimal 
diversity sum 0.8660. The second constellation is 3 dimensional constellation which has 5 
elements {R = 0.7740) with diversity product 0.7183 and diversity sum 0.7454. The third 
constellation is a 4 dimensional one consisting of 9 elements {R = 0.7925) with diversity 
product 0.5904 and diversity sum 0.6403. Here we used Genetic Algorithm to optimize the 
diversity product and the diversity sum at the same time to acquire the last two constella¬ 
tions. 

One can see that around 5 db, the second constellation surpasses the hrst one and is 
getting better and better as SNR becomes larger. This can be easily understood that the 
diversity function of the hrst constellation is approximately dominated by 1/p^ at high SNR, 
while the diversity function of the second constellation is dominated by 1/p®. The same 
explanation can be applied to the third constellation’s performance. One can even foresee 
that higher dimensional constellations will perform even better and the BER curve will be 
sharper than low dimensional ones. It is believable that higher dimensional constellation 
ones will have much more diversity gain compared to lower dimensional ones as for other 
transmission rate. 
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Figure 5: Performance of different dimensional constellations with the same rate 

4.4 Constellations optimized for certain values of SNR by diver¬ 
sity function 

Different industrial applications require different level of reliability of the communication 
channels. One may want to optimize the constellation at certain Block Error Rate (BER) 
or Signal Noise Ratio (SNR). Algebraically designing codes for this purpose seems to be too 
complicated, while the numerical method can be easily applied for this purpose. In fact as we 
mentioned before, algebraic structure is still benehcial for optimizing the diversity function 
using simulated annealing algorithm, we also hnd our genetic algorithm only needs minor 
modification to £t in for this purpose. We have the following graph to show some further 
constellations. 

We have already seen the first two constellations, the third constellation is specifically 
optimized at 25 db using Simulated Annealing Algorithm based on structure , one can 
see that optimizing at different curve does change the shape of curve. Consequently one can 
see the flexibility of the numerical method employed here. 

The third constellation has the following representation: 


A^B^\A 


-0.3018715 +0.8863567 0.1337423 - 0.3245896 * i \ 

-0.3487261 - 4.0441897E - 02 * i -0.9215508 - 0.1658271 *i ) ' 


B 


-0.9599144 - 0.1577551 * i 0.2074585 + 0.1031435 * i 
-0.2074221 + 0.1032166 * i -0.9598588 + 0.1580936 * i 
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Figure 6: Optimize the constellations at certain SNR by diversity function 

4.5 General form constellation numerical design 

As above, we discussed numerical methods for the special form of unitary constellation: 


where 


a. ^ 




There are numerical methods available to design general constellation as well. Right now 
we are only taking the simplest approach. In this case, we assume that T = 2M and we 
consider the following constellation: 


{A'^B\A e U{M),B 





With this algebraic structure, whenever we want to calculate the diversity product or sum 
or more generally the diversity function, one needs to calculate It is easily verihed 

that they will end up being one of the following: 


B*A^B k = 

In this way, we have a hnitely generated constellation with reduced number of targets. 
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Table 10. The following tables show our elementary results using SA. 


size 

3 

4 

5 

6 

7 

8 

9 

rate 

0.3962 

0.5000 

0.5905 

0.6462 

0.7018 

0.7500 

0.7925 

diversity sum 

0.8654 

0.7901 

0.7889 

0.7652 

0.7514 

0.7422 

0.7369 


size 

3 

4 

5 

6 

7 

8 

9 

rate 

0.3962 

0.5000 

0.5905 

0.6462 

0.7018 

0.7500 

0.7925 

diversity product 

0.8582 

0.7424 

0.7330 

0.6450 

0.6361 

0.6216 

0.5822 


5 Conclusion and Future Work 


In this paper, we study the limiting behavior of the diversity function by either letting the 
SNR go to inhnity or to zero. Several criteria for designing unitary constellations are derived 
from the analysis of the limiting behavior. Two numerical algorithms, SA and GA produce 
excellent codes with large diversity product and diversity sum. The numerical methods can 
equally be employed to optimize the diversity function at a certain SNR. 

Some algebraic structured constellations are presented. Is there any other or better 
constellation structure? In |^, parametric codes are designed using a structure. 


Can we also design other codes algebraically using other structure? It should be an interesting 
question worth exploring further. 

One can see that the numerical algorithms produce excellent codes frequently although 
only some simple random mechanism is applied. To explain this, we need explore further the 
deeper geometric structure of U{M) and the Stiefel manifold St,m- Also our implementation 
of the algorithm is very simple and naive. We are still wondering if more subtle and delicate 
implementation will generate more remarkable results. 

We already have some examples showing that the diversity product (or sum) doesn’t 
necessarily mean good performance. Basically they are just the maximized minimum dis¬ 
tance. We know that some of the turbo or LDPC codes don’t have large minimum distance 
either, but they are capacity achieving codes which perform very well. The reason for this is 
because these codes have remarkable distance distribution. Is there any possibility to look 
further into the transformation of the criterion? For instance, consider the sum of all the 
distance. In this case, better distance distribution seems to be guaranteed. In this case the 
problem can be converted to a polynomial minimization problem which is widely explored 
in optimization theory. 

We did some elementary explorations for general type constellations. Compared to the 
special form of constellations, general form constellations don’t have the unitary differential 
modulation which will speed up the transmission rate. Maybe it is possible to develop a 
differential modulation scheme for general constellations? What performance could a general 
form constellation have compared to a special form one? It seems like that a lot of further 
research will be needed to address these questions. 
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Appendix 


A An Upper Bound for the Optimal Diversity Product 
and Sum for 3 Element Constellations 


In this appendix, we are going to give an npper bound for the diversity product of an arbitrary 
dimensional constellation with exactly three elements. As a consequence, one can have the 
optimal value and optimal conditions for a 2 dimensional constellation with 3 elements. The 
same technique can be used to give the optimal bound for diversity sum. The main idea 
follows from |^. 

Let A = (aij)nxn he a n x n matrix and denote with Sn the symmetric group in n 
elements. Let 

n 


F{n) 


sup 

AeU(n} 


E n 

uGSn ^=1 




Then one has: 


Theorem A.l. Let V be a unitary constellation as a subset ofU{n) with 3 elements, then 

sup JJV < (A.l) 

VcU{n) ^ 

For the proof of this theorem we will need the following technical lemma: Let S{m,n) 

n 

be the set of n x m matrices of the form d) = where 0 < < vr and ^ = tt. Let 

i=l 

m 

di{^) := P sin^jj, then we have: 

Lemma A.2. 

/ 7r\ ™ 

max min = sin — ) 

<I>eS(m,n) ,n V n/ 

and equality holds if and only if ^ij = ^, i = 1,2, ■ ■ ■ ,n and j = 1,2, ■ ■ ■ ,m. 

Proof. Let d be the left hand side in above formula. Then 


d” < n ^ =n n ® ^=n n 

2=1 2=1 j = l j = l 2=1 

E*=i sin ^ij 


1=1 


n 


1=1 


that means 


, , 7r\™ 

a < I sm — ) . 

nJ 


On the other hand if one chooses = - one sees that 

n 

( ^ 
d > sm — 

V n 

which establishes the claimed equality. We leave it to the reader to verify that the maximal 
value is achieved in a unique manner. □ 
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Proof of Theorem \A.1\ . Consider a constellation V C U{n) with cardinality |V| = 3. Withont 
loss of generality we can assnme that V has the form: 


where is the n x n identity matrix and D is a diagonal matrix of the form 

D = diag 

and A is an arbitrary nnitary matrix. Assnme A has eigenvalnes ■ ■ ■ > i-e. there is 
a nnitary matrix U = {uij)nxn such that 

U-^AU = diag ■ • • , 

If either 

I det(/ — D)I < 3^ or I det(/ — A)| <3^ 
then automatically we have 

Assume therefore that | det(/ — D)\ >3^ and | det(/ — A)| > 3t, that is 

|det(J-diag(e'''^e'"^••• ,e'®"))| > 3?, 

|det(/-diag(e*^^e*^^•■■ ,e*‘^")| > 3?, 
then according to Lemma |A.2| we have the following inequality: 


det(D-A)| = |det(D-f/diag(e*‘^Se*‘^=,--- , 

= I det{DU - t/diag • • • , e*‘^"))| = | det{{uij{P^^ - P^^)))\ 

n 

< E n 

cr&Sn i=l 


Taking the nth root and dividing the result by 2 shows also in this case that the diversity 
product is at most the value on the right hand side in (|A.1|) . □ 

In general we do not know how sharp the estimate is. However when n = 2 we have: 

Theorem A.3. Let V be a unitary constellation as a subset ofU(2) with 3 elements, then 


sup TT V = a/3/2. 

vcn(n) 


Proof Since U{2) has the explicit form (|4.1|) one can verify that F{2) 
rem A-l| , we have 


1. From Theo- 


sup TT V < v^/2. 

vec2,|v|=3'‘"'‘" 


On the other hand Remark A-4| shows the existence constellations of three elements whose 
diversity product is a/3/2. □ 
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Remark A.4. If V C U{2) is a constellation of cardinality 3 having then V 

must have the following form: 


{C,CADA-\CBEB-^} or {C,CAFA-\CBGB-^} or 
{C,ADA-^C,BEB-^C} or {C, AEA'^C, BGB-^C} 

where 



and A, B,G are arbitrary 2x2 unitary matrices. 

Remark A.5. If n = 3 then we computed with the help of Matlab E{3) = 1.299. This 
results in the upper bound: 


sup TT V < ^E{3)V3/2 = 0.95 

VCf/(n) 

When n > 4, we believe E{n) > (:^)"' and the inequality should be trivial. 

For the diversity sum of a 2 x 2 constellation with 3 elements, we have the same result 
as for the diversity product. The following Lemma deals with diagonal unitary constellation 
hrst. 


Lemma A. 6 . Let V be a diagonal unitary constellation as a subset ofU{2) with 3 elements, 
then we have 

sup ^2 ^ = v^/2, 

VcD{n) 

where D{n) denote the set of all the diagonal unitary matrices. 

Proof. Without loss of generality, we assume that V = {I 2 , A, B}, where /2 is 2 x 2 matrix 
and 



Suppose by contradiction 

sup ^ V > V^/2, 

VcD{n) 

that is there are matrices A, B such that. 



- e*®" 

2 

> 

6, 

\l-e^^^f + \1 


2 

> 

6, 


- e*®" 

2 

> 

6. 
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One can verify that after changing of variables, the three above equation can be rewritten 
in the following way: 

sin^ X + sin^ a:i > ^, sin^ y + sin^ l/i > ^ and sin^ z + sin^ Zi > ^, 

where ^ < x,y, z,Xi,yi, Zi < x + r/ + z = tt and Xi + yi + Zi = tt. 

By adding up three inequalities one can derive: 

sin^ X + sin^ y + sin^ ^ + sin^ Xi + sin^ yi + sin^ ^ 2 (^•^) 

Since the function sin^(x) is concave on the interval [|, |] one shows that 


sin^ X + sin^ y + sin^ 2; + sin^ Xi + sin^ yi + sin^ Zi 


.0 • 2^ + y + z , Q • 2 2^1 +2/1 + 0 

< o sin-h 6 sin - 


which contradicts ( A.2 ). 

On the other hand, let 

u = ll,A,B\A = 
one verihes that 


0 

- 2 ni \ , n — 
0 6 3 


6 3 0 

-47ri 

0 6 3 


= V 3 / 2 . 

That means the upper bound can be reached. 


9 

2 ’ 


□ 


With the lemma above, we can prove the following theorem. 

Theorem A.7. Let V he a unitary constellation as a subset ofU{2) with 3 elements, then 
we have 

sup y'V = v/3/2. (A.3) 

VzU(n) 

Proof. Consider a constellation V C f/(2) with cardinality |V| = 3. Without loss of generality 
we can assume that V has the form: 


V = {J2,Z^,A}, 

where is the n x n identity matrix and D is a diagonal matrix of the form 

D = diag (6*+, 6*®^) 

and A is an arbitrary unitary matrix. Assume A has eigenvalues 6*‘^b6*‘^2, i.e. there is a 
unitary matrix U = {uij) 2 x 2 such that 

U-^AU = diag ( 6 + 1 , 6 +=). 
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If either 

then automatically we have 


D\\^p < 6 or \\I - A)\\p < 6 

V3 


Ev 


< 


Assume therefore that \\I — D 11^ > 6 and ||/ — y4)||^ > 6, that is 

||J-diag(e*^\e*®=)||?,>6, 

||/-diag(e*‘"\e*^^)||^>6, 

then according to Lemma A.6 and the explicit form ^ of U{2), we have the following 
inequality: 


\D-A\\l = \\D-Udiag{e^‘^\e^^^)U-^\\l = ||Df/- f/diag 


= \\[Uij(e 


idi 


^vpj' 


I | 2 n 

= |mii| (|e 


+ + M") = 6. 


idi _ I 

^^2 


+ e 


*02 


3*¥52 


') 


Taking the 2th root and dividing the result by 2^2 shows also in this case that the diversity 
sum is at most the value on the right hand side in ( |A.3 ). 

On the other hand one can take the same constellation as in Remark |A.4| and show that 
the upper bound for the diversity sum is reached. □ 
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